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3 . A method of halftoning according to claim 1 , wherein: 

the next scanline error impulse response is a sampling of a Cauchy distribution, 
said sampling being normalised so that a sum of next scanline error impulse response 
values is unity. 

5 

4. A method of halftoning according to claim 1, wherein: 

the next scanline error impulse response is left-right symmetric. 

5. A method of generating an error diffusion mask suitable for use with any of the 
10 aforementioned methods. 

6. An error diffusion mask suitable for use with any of the aforementioned 
methods. 

15 7. A method of halftoning an image, said image comprising a plurality of pixels 
each having an input value and an assignable output value that can take on one of at least 
two output values, where pixels are processed scanline by scanline, and scanlines are 
processed one at a time from the top of the image to the bottom of the image, and where a 
scanline is processed pixel by pixel either from left to right or from right to left, and 

20 where the processing for each pixel comprises the steps of: 

(a) determining the output value of a current pixel using a sum of the input 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 

25 for the pixel, and secondly the output value of the pixel; 
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(c) adding proportions of the error at the current pixel to the neighbourhood 
error values of yet to be processed pixels of the current and next scanline; 

and where the said proportions of the error at a current pixel are designed so that 
the next scanline error impulse response, being that function which maps 
5 (A) from a horizontal pixel offset; 

(B) to the total proportion of the error at the current pixel added to the 
neighbourhood error of that pixel of the next scanline which is displaced by the horizontal 
pixel offset from the current pixel, following complete processing of the current scanline; 

approximates a function which spreads with self-convolution in proportion to the 
10 degree of self-convolution. 

8. A method as claimed in claim 7, where the next scanline error impulse response 
approximates a scaled sampling of a Cauchy distribution. 

15 9. A method as claimed in claim 7, where the next scanline error impulse response 
approximates a function which has a Discrete Space Fourier Transform which is a 
replicated two-sided exponential function. 

10. A method as claimed in claims 7 to 9, where the output value of a current pixel is 
20 determined by comparison of the sum of the input value of the current pixel and the 

neighbourhood error value for the pixel against a threshold value. 

11. A method as claimed in claim 10, where in step (c), for a current pixel at pixel 
position (ij), being column i and scanline j, error at the current pixel is added to the 

25 neighbourhood error of only those pixels which are either: 
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(i) on the current scanline ahead of the current pixel at a pixel position 
(i+current__offset, j), where, for left to right processing of the current scanline, 
current_offset is greater than zero, and, for right to left processing of the current scanline, 
current_offset is less than zero, or 
5 (ii) on the next scanline below or behind the current pixel at a pixel 

position, (i-next_offset, j+1), where, for left to right processing of the current scanline, 
next_offset is greater than or equal to zero, and, for right to left processing of the current 
scanline, next_offset is less than or equal than zero. 

10 12. A method of halftoning an image, said image comprising a plurality of pixels 
each having an input value and an assignable output value that can take on one of at least 
two output values, where pixels are processed scanline by scanline and scanlines are 
processed one at a time from the top of the image to the bottom of the image, and where a 
scanline is processed pixel by pixel either from left to right or from right to left, and 

15 where the processing for each pixel comprises the steps of: 

(a) determining the output value of a current pixel using a sum of the input 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 

20 for the pixel, and secondly the output value of the pixel; 

(c) selecting, using the current pixel input value, a set of proportions and a 
set of corresponding pixel position offsets, from a family of sets of proportions and 
corresponding pixel position offsets; 
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(d) adding the selected proportions of the error at the current pixel to the 
neighbourhood error values of yet to be processed pixels at pixel positions offset from the 
current pixel by the selected corresponding pixel position offsets; 

and where each set of the said family of sets of proportions and corresponding 
5 pixel offsets, is designed so that the next scanline error impulse response corresponding to 
that set, being that function which maps 

(A) from a horizontal pixel offset; 

(B) to the proportion of the error at the current pixel added to the 
neighbourhood error of that pixel of the next scanline displaced by the horizontal pixel 

10 offset from the current pixel, following complete halftone processing of the current 
scanline using only the said set of proportions and corresponding pixel offsets; 

approximates a function which spreads with self-convolution in proportion to the 
degree of self-convolution. 

15 13. A method as claimed in claim 12, where each next scanline error impulse 
response, corresponding to a set of proportions and pixel offsets, approximates a scaled 
sampling of a Cauchy distribution. 

14. A method as claimed in claim 12, where each next scanline error impulse 
20 response, corresponding to a set of proportions and pixel offsets, approximates a function 

which has a Discrete Space Fourier Transform which is a replicated two-sided 
exponential function. 

15. A method as claimed in claim 13 or 14, where the family of sets of proportions 
25 and pixel offsets together with the selection, using the current pixel input value, of a set of 
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proportions and pixel offsets, are designed so as to minimise processing while also 
minimising the presence of artifacts in the halftone output artifacts including cross-over 
artifacts and poor spreading in sparse halftone patterns. 

5 16. A method as claimed in claim 15, where the maximum absolute offset in each set 
of the family of sets of proportions and pixel offsets, varies so that the family of sets 
includes a set with small maximum absolute offset and a set with large maximum ■ 
absolute offset, and where intermediate input values primarily select sets with small 
maximum absolute offset, and extreme input values primarily select sets with large 
10 maximum absolute offset. 

17. A method as claimed in claim 16, where the output value of a current pixel is 
determined by comparison of the sum of the input value of the current pixel and the 
neighbourhood error value for the pixel against a threshold value. 

15 

18. A method as claimed in claim 17, where in step (d), for a current pixel at pixel 
position (ij), being column i and scanline j, error at the current pixel is added to the 
neighbourhood error of only those pixels which are either: 

(i) on the current scanline ahead of the current pixel at a pixel position 
20 (i+current_offset, j), where, for left to right processing of the current scanline, 

current_offset is greater than zero, and, for right to left processing of the current scanline, 
current_offset is less than zero, or 

(ii) on the next scanline below or behind the current pixel at a pixel 
position, (i-next_offset, j+1), where, for left to right processing of the current scanline, 
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next_offset is greater than or equal to zero, and, for right to left processing of the current 
scanline, next_offset is less than or equal than zero. 

19. A method of halftoning an image, said image comprising a plurality of pixels 
each having an input value and an assignable output value that can take on one of at least 
two output values, where pixels are processed scanline by scanline and scanlines are 
processed one at a time from the top of the image to the bottom of the image, and where a 
scanline is processed pixel by pixel either from left to right or from right to left, and 
where the processing for each pixel comprises the steps of: 

(a) determining the output value of a current pixel using a sum of the input 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 
for the pixel, and secondly the output value of the pixel; 

(c) selecting, using the current pixel input value, a set of proportions and a 
set of corresponding pixel position offsets, from a family of sets of proportions and 
corresponding pixel position offsets, said set of proportions being in accordance with a 
next scanline error impulse response that approximates a function which spreads with 
self-convolution in proportion to a degree of self-convolution; 

(d) adding the selected proportions of the error at the current pixel to the 
neighbourhood error values of yet to be processed pixels at pixel positions offset from the 
current pixel by the selected corresponding pixel position offsets, 

and where each set of the said family of sets of proportions and corresponding 
pixel offsets, only includes pixel offsets corresponding to pixels on the same scanline as 
the current pixel or to pixels on the next scanline. 
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20. An apparatus for halftoning an image, said apparatus comprising: 

means for determining an output value of a current pixel on a current scanline 
using a sum of an input value for the current pixel and a neighbourhood error value at the 
current pixel; 

means for determining an error at the current pixel as the difference between (i) 
the sum of the input value for the current pixel and the neighbourhood error value at the 
current pixel, and (ii) the output value of the current pixel; and 

means for adding a proportion of the error at the current pixel to neighbourhood 
error values at as yet unprocessed pixels of a subsequent scanline in accordance with a 
next scanline error impulse response; wherein said next scanline error impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 

21. An apparatus for halftoning an image, said apparatus comprising: 
a memory for storing a program; 

a processor for executing the program, said program comprising: 

code for a determining step for determining an output value of a current pixel on 

a current scanline using a sum of an input value for the current pixel and a neighbourhood 

error value at the current pixel; 

code for a determining step for determining an error at the current pixel as the 

difference between (i) the sum of the input value for the current pixel and the 

neighbourhood error value at the current pixel, and (ii) the output value of the current 

pixel; and 
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code for an adding step for adding a proportion of the error at the current pixel to 
neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
accordance with a next scanline error impulse response; wherein said next scanline error 
impulse response: 

5 approximates a function which spreads with self-convolution in proportion to a 

degree of self-convolution. 

22. A computer program product including a computer readable medium having 

recorded thereon a computer program for directing a processor to execute a method for 
10 halftoning an image, said program comprising: 

code for a determining step for determining an output value of a current pixel on 

a current scanline using a sum of an input value for the current pixel and a neighbourhood 

error value at the current pixel; 

code for a determining step for determining an error at the current pixel as the 
15 difference between (i) the sum of the input value for the current pixel and the 

neighbourhood error value at the current pixel, and (ii) the output value of the current 

pixel; and 

code for an adding step for adding a proportion of the error at the current pixel to 
neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
20 accordance with a next scanline error impulse response; wherein said next scanline error 
impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 
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23. A computer program for directing a processor to execute a method for halftoning 
an image, said program comprising: 

code for a determining step for determining an output value of a current pixel on 
a current scanline using a sum of an input value for the current pixel and a neighbourhood 
5 error value at the current pixel; 

code for a determining step for deteimining an error at the current pixel as the 
difference between (i) the sum of the input value for the current pixel and the 
neighbourhood error value at the current pixel, and (ii) the output value of the current 
pixel; and 

10 code f or ^ adding step for adding a proportion of the error at the current pixel to 

neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
accordance with a next scanline error impulse response; wherein said next scanline error 
impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
15 degree of self-convolution. 
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Patent Attorneys for the Applicant 
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